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ABSTRACT 


A ccuputer model was develored to simulate a 
Frojectile ram-pfropellant feed system for a lLiguid 
Propellant Gun. Using a lumped parameter approach, a 
set of sinultaneous differential equations was derived 
for the ccmplex interaction of the propellant fluid, 
the drivirg injector and the projectile. The computer 
model was verified against a 20 mm experimental 
apparatus. Injector displacement, projectile 
displacement, and chamber pressure were ccmpared for a 
nominal driving pressure of 140 psi. The important 
system farameters affecting projectile ram time and 
chamber fressure cscillations were investigated and 
potential prcktlem areas for testing with actual 
Fropellant were identified. 
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I. INTRODUCTION 


As the state cf the art of weapons technology improves, 
mcre complex systems are continually being developed. To 
neet future needs these systems must have improved 
performance, less weight and volume, less conplicated 
logistics, and at the same time be more cost effective than 
current systems. One candidate for these future weapon 
systems is the liguid propellant gun (LPG). Nuperous 
investigaticrs over the past thirty years have demonstrated 
its potential if the technology can be developed to utilize 
a liguid frogellant in a large caliber gun. 


One impertant area that needs investigation is the 
aspect of fluid prcpellant handling. The ordnance designer 
must kncw the necessary propellant flow rates and pressure 
regimes needed to achieve the desired system performance. 
These parameters must be within allowable safety margins for 
the handling and tse of the explosive fropellant. One 
possibie design of an LPG which would allow a very high 
firing rate utilizes the propellant to ram the loaded 
projectile. Ccmbining the normally separate fropellant load 
and ram cycles decreases the required sovement of mechanical 
parts. Thus, using this method, it is feasible that firing 
rates in excess of four to five times present rates could be 
achieved. This thesis was directed toward identifying the 
impcrtant flvid dynamic parameters involved with such a 
prcjectile ram feed systen. 





II. BACKGROUND 


—_ oo or oe oe ee oe oe ee 


During 16°74 and early 1975, an investigation was 
conducted at the Naval Postgraduate Schocl to study, both 
analytically and experimentally, the fluid dynamics of a 
liguid propellant under conditions Similar to those which 
would exist in a rapid-fire LPG feed system. The results of 
this investigation were to be used to establish such LPG 
design and performance parameters as time-to-load, injection 
SuUfFly fressure, injection system configuration, ullage, 
charge-tc-mass ratio, caliber size, and projectile mass. It 
was hoped that this investigation would identify any 
Fetentia] preblem areas for further detailed research. 


The kasic objective of the experimental portion cf the 
investigaticn was to identify what fluid dynamic 
characteristics of a liquid fpropellant feed system would 
limit lcading times and hence rates of fire. A 2zO aon 
experimental model of a kasic propellant feed system was 
designed and ruilt. Data on injector displacement, treech 
chamber fressure, and ram gas pressure (input driving 
Fressure), were reccrded for driving fressures between 50 
and 220 psig in 10 psi increments. As reported in Ref. 1, 
it was fcund that the instantaneous behavior of the chamber 
pressure was the result of a complex interaction of inertia 
ferces, visccus forces, and the unsteady motion of the 
fluid. Ike experiments demonstrated that frictional and 
inertial effects were significant during the movements of 
the injectcr and the projectile slug. Once the projectile 
slug stopped, the effect of entrapped gas in the fluid 
caused large breech chamber pressure oscillations. On 
several runs, sSub-atmospheric pressures were experienced 
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which suggested the possibility of cavitation and hence 
vapor-phase igniticn. Over all, it was found that the ran 


tize had a quadratic dependence upon ram pressure. 


The analytical portion of the investigation was directed 
toward fredicting the pressure and flcw rate of the 
prepellant ard the projectile slug motion during an LPG 
loading cycle. The mathematical model which was developed 
and programmed on @n analog computer predicted the pfesition, 
velocity, and acceleration of the projectile anc the 
Fressure at various points in the system as functions of 
time. The mcdel was tested against experimental results and 
found to ke adequate for the prediciton of projectile can 
time. An analysis was performed using the model to indicate 
the areas cf system redesign likely to be most prefitable 
and to cktain preliminary predictions of LEG loading systen 
performance under a variety of design ccnditions. The 
results cf tke analytical study, as well as a summary cf the 


experimental work is contained in Ref. 2. 


The afcrementioned model was the starting point for the 
present study. To better understand the derivations in the 
following secticns, Figure (1) indicates the geometry cf a 
kasic prce-ectile ram fropellant feed system. An injector 
chamber is filled with propellant, which is then pumped into 
the gun -reech by applying ram pressure to the ram side of 
the injectcr piston. This can be accomplished by using high 
fressure cas from an accumulator or by using an hydraulic 
drive systen. The force exerted by the injector fiston 
drives the fluid propellant through the connecting line into 
the gun kEreech. The rapidly accelerated propellant drives 
the already lcaded projectile to its seated positicn ready 
for firing. 
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Figure 1 - BASIC GEOMETRY 





The ccuplex nature of this dynamic ram process is not 
easy to describe with mathematical equations. Many 
researchers have approached the froblem using wave mechanics 
and partial differential eguations. These eguations were 
then modified according to the propagation characteristics 
of the system and the various boundary conditions 
encountered therein. This usually led to involved finite 
difference methcds of soluticn on a digital computer. 


Another fossikle approach considers the kinetic energy 
invclved in the ram frocess. This is the apfroach that was 
used to derive the governing eguations for the previously 
described analcg computer model. Beginning with the input 
side, a force balance on the injector piston yielded an 
eguaticn fcr the pressure in the injector chamber in terms 
of the injector motion. Writing Bernoulli's equations for 
head losses at the expansion to the breech chamter, the 
ccntracticn tc the connecting line, and an orifice in the 
connecting line resulted in an equation fcr the fressure 
drcp in tte connecting line in terms of the square cf the 
fluid velccity and the fluid acceleration. A force balance 
on the freiectile slug yielded a third eguation which 
described tke Ereech chamber pressure in terms of the slug 
motion. Ic solve these equations, scme method of relating 
the moticns in the injector, connecting line, and _ tkreéch 
chamber had to tke established. It was assumed that the 
fluid was inccmpressible and therefore these moticns were 
identical. Ic faciditate the analog simulation it was also 
assumed that the irput ram pressure was a step input. In 
the process of converting these eguations to an integral 
form for wiring on an analog computer, all variables were 


ncermalized tc permit scaling. 


Despite these limiting assumptions, the analog cogputer 
mcdel was able to predict the time to ram the projectile and 
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the chamrer ;ressures which cccur during the time that the 
injectcr is in motion. Due to the incompressitility 
assumpticn, the model was unable to account fcr the pressure 
transients cr "water hammer" type pressure oscillations 
which occurred when the projectile was seated and the 
propellant fluid was being decelerated. These pressure 
oscillaticns are important to the system designer because 
the peak kEreech chamber load pressure is experienced during 
these oscillations. It is also possible that these 
oscillaticns could interfere with the uniformity cf the 
ignition and subsequent combustion of the propellant. 


A straight fcrward attempt was made to extend this 
analog mcdel Ty first converting it to a digital computer 
prcgranm. This allowed an accurate representation cf the 
input ram pressure and the use of unnormalized variables to 
be incorrperated. These improvements increased the accuracy 
of the mcdel tut the effect of the compressibility of the 
prepellant fluid still was not taken into account. Irying 
to add compressibility effects to this model by adding time 
derivatives cf the pressure terms to the gcverning equations 
unlinked tte motions of the injector, the connecting line 
fluid, and the frojectile slug causing a problem with too 
Many unkncwn variables for the number of equations invclved. 


It was decided that a new approach should be tried which 
would adeguately describe the system pressure oscillaticns. 
The apprceach would feature an engineering mcdel which would 
be easy to use, adaptable to any LPG feed system, and not 
obscure tke interaction of system parameters by complicated 
mathematics. 
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Tit. DERIVATION C 


A. LUMFEC EXRRAMETER APPROACH 


The apfrreach that showed the greatest promise for 
modeling the LPG feed system was that of the pati alys| 
transmission line cencept. This approach, which has kLecome 
FCfular in tke last five to ten years, is based on a 
pressure-vcltage and flow velccity-current analogy with 
Electrical Engineering determinations. It is an outcrowth 
ef the larce amcunt of effort that has been devoted to 
investigating fluid line transients. Reference 3 is a_ good 
survey of this field. 


For ccmplex systems this approach is usually simplified 
ky using tte approximation of lumped farameters. The 
effects cf fluid inertia, capacitance, and resistance are 
"lumped" and considered to act only in discreet areas cf the 
system. This results in a reduction cf the unknown 
parameters due tc the assumed lack of interaction cf the 
different fluid effects. This approximation, however, 
results in the mecessity of using several empirical 
ccnstants which must be determined by fitting experinental 
Gata. The crdinary differential equations that are derived 
frcm this method can be solved either by Laplace 
transfcrmaticn or Ly computer integration... 


To censider the effect of fluid inertia it is agsumed 


that cnly fressure and inertia forces are present anc that 


ccmpressibility effects are negligible in the volume under 
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consideraticr. Then, 
dV 
P -P = pLlL— 
12 oe 
where 0 is the propellant Bidaiad density, -L ais “a 


characteristic length and ae because the flcw is 


considered irccmpressible for this building block (Fig.2A). 


To ccnsider the effect of fluid capacitance, it is 
assumed that only compressibility effects are important, and 
that inertia and resistance effects may be neglected in the 


vclure under consideration. Therefore, 
V V Evo 
1 2 Reedy 


where K is the effective bulk modulus cf the propellant 
fluid, L is a characteristic length and ett (Fig. ZB). 


Eecause fluid resistance can be affected by sc many 
different parameters, it is impossible to write a general 
equation describing the pressure drop due to fluid 
resistance. It is best to treat it empirically using an 
experimentally derived figure. Thus, in the volume under 


consideraticr, 


P -P =RVYV 
] 2 ee. 

where Rv is a function of fluid velocity and ee (Fig. 
2C). Changes in cross sectional area can be accounted for 
by using apprcfriate area ratios. 

These results specify three building biccks which can be 
ccubined in any sequence to nodel the dynamic 
characteristics of a system. The complexity of the model 


can ke increased to any degree necessary by including more 
and more ccmkinaticns of these three basic building blocks. 
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A) INERTIA 


| dV — 
Lee 7 one 
B) CAPACITANCE 
V V 
1 2 
1 2 
ap a Pp =P =P 
en: Fy | ole 
C) RESISTANCE ~ 
1 2 
P. P, 
P -~P =Rv V VY =\Ve= vy 
1 2 1 2 
Figure 2-—~- LUMEED EARAMETER MOLULES 





Be. LPG SYSTEM GOVERNING EQUATIONS 


Several aifferent combinaticns of buildirg blocks’ were 
tried in attempting to model the experimental LPG systen. 
The experimental results showed that inertial and resistance 
blocks should Ee included in the model. The fact that 
Fressure transients occurred in the breech chamber indicated 
that a capacitance block should also te included. To 
ninimize tte use of computer time, it was decided to start 
with the mcst simple model thought adequate and work toward 
more complex models as the quality of the computer results 
dictated. 


The first attempt considered that ccmpressibilty effects 
would dcmrinate in the injector chamber as the accelerating 
piston interacted with the propellant fluid. The resistance 
effects cf changes in the feedline cross sectional diameter 
and associated orifices were modeled as one resistance block 
at the keginning of the connecting line. Inertial effects 
were thought to dcminate in the connecting line. This would 
be particularly true in large scale models where the length 
of the ccnnecting line would be very large. It was thought 
that compressibility effects should again dominate in the 
kreech chamber as the accelerating propellant fluid drove 
the projectile dcwn the breech chamber. Results in this 
case were ocverned by inertial effects because the 
calculated ccnnecting line velocity was too great. This led 
to a decrease in the chamber pressure as the projectile 
accelerated until a large negative chamber fressure existed. 
In additicn to this deficiency, the computed system fressure 


cscillaticns were not sufficiently damped. 


To try tc adjust the connecting line velocity, ancther 
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resistance rlilcck was added. As will be seen in the section 
discussing the parameters which affect the pressure 
oscillaticns, the value of this resistance coefficient is 
prcrorticral to the amount of damping in the fressure 
oscillaticns. Even with a value which gave the froper 
damping, tke breech chamber pressure was dominated Ly the 
Fressure loss due to the accelerating fluid in the 
ccnnecting line. 


To reduce the dominance of the fluid inertia effects, it 
was decided tc restructure the model. fo mcdel the intector 
chamber it was decided to consider the effect of fluid 
capacitance anc inertia. The connecting line was modeled by 
twce resistance blocks separated by a fluid capacitance 
bicck. The kreech chamber was modeled using a fluid 
capacitance and a fluid inertia block, as was done for the 
injectcr chamber. This mcdel gave the final results 
described in the section comparing computer and experinental 
cata. Ey ccmbining these building biocks, the governing 
egtaticns became: (See Fig. (1) for notaticn) 


For the intectcr: 


1, PRA -PA =MpV. +Kes V+ Por A 
1 ie 1 le 1 


ae ep 
ot 2 em 
Sc al ey ll 


For the ccnnecting line: 
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These equations were subseguently mcdified tc acccunt 
for changes in cross-sectional area and changing gecmetry as 
the injectcr chasber decreased in volume and tne kreech 
chamber increased. Once the projectile reached the end of 
the breech ckamber, the force balance equaticn was nce lcnger 
considered applicable. At this time, it was noted that the 
remaining ecuations for the connecting iine and the rreech 
chamber cculd be combined to forn a second crder 
differential equation for the breech chamber pressure which 
descrikred tre "water hammer" type pressure oscillations 


displayed by the experimental data. Thus, 





ut A3 RV> p % K K 


10. += —— P 
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Once tke injector piston reached the end of its travel, 
the injectcr chamber no longer existed and the gcverning 
egtaticns were considered as no lcnger applicatkle. 
Ultimately, the injector chamber pressure was prcepagated 
down the ccnnecting line as the propellant fluid came to 
rest. Thrcuchout these analyses, input ram pressure was 
modeled using expcrential terms to fit experimental data. 


The eguaticns frem the lumred parameter approach fcr the 
Frepellant fluid were combined with the force balance 
eguations fer the injector piston and the frojectile using 
State Variatle methods. Subsequently, a computer program 
was constructed to solve these state variable relations 
using a forrth crder Runga-Kutta integration routine for 
sinultanecus first order differential equations which was 
developed at the Naval Postgraduate School. A brief 
descripticn cf the state variable method and a listing of 
the ccmputer program can be found in Appendix A. 


The fcllcwing is a listing cf the nomenclature and the 
vaiuves cf ccnstants used in the derivation cf the governing 
equations anc the computer program: 


2 2 
A1 Injectcr 1.77 in 11.42 cr 
Crose-Sectional Area 
; : 2 2 
Az Ccnnectirg Line «255 rn 1.65 c# 
Cross-Sectional Area 
ee 2 
A3 Ereech Chamber 49 in 3.16 cB 
Cross-Sectional Area 
“eZ 2 
AR Ram Fistcn 1.77 in 11.42 co 
Crose-Sectionéel Area 
11 Injectcr Length 1.60 in 4.22 cr 


So a ee me SS ae a eras mia 





£2 Connecting Line Length 16.0 in 40.6 cr 


L3 Breech Chamber Length 5.0 in - 12,7 <¢e 
ME Injectcr Piston Mass -0052 meee 913.8 gt 
in 
MS Projectile Mass -00053 wwsee 93.1 gs 
in 
KFP Injectcr Friction -001 ib-sec 175.7 gny/sec 
Factcr in 
KFS Frojectile Friction 001 ib-sec 175.7 guysec 
Factcr in 
EDI Injector Back Pressure 24.0 PSI 163.3 kEa 
FOS Projectile Back -50 FSI 3.4 kEa 
Pressure 
K Fffective Bulk Mcdulus 3200 PSI 6808 kEa 
RHC Propellant Density -000093 apes 994 gee 
_ 4 
in 
PK Ram Fresesure Ccmputed 
P1 Injectcr Chamber Pressure Conputed 
v1 Injectcr Pistcn Velocity Ccmputed 
V2 Injector Chamber Exit Velocity Ccmnputed 
PY Connecting Line Pressure Ccmputed 
V6 Breech Chamber Entrance Velocity Ccmputed 
F8 Breech Chamber Pressure Computed 
v8 Projectile Velocity Computed 


The treatment cf several areas of the computer mnecdel 
continually reappeared as reguiring refinement. The 
governing differential equations worked well during the 
dynamic ferticn cf the feed cycle but experienced 
difficulties during the initial and final static periods. 
The handling cf static friction and back pressure as a 
ccrstant valve created the possibility of negative 
velocities urtil the driving ram pressure overcame the 
systen kack pressure. These negative values never occurred 
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in the real system due to the gecmetric restraints on the 
injector fiston and projectile; therefore, the computer 
model had tc be manipulated to maintain this ccndition. 
Unfortunately the need for simultaneous soluticn cf the 
governing equations made it difficult to manipulate the 
initial ccnditions without greatly increasing the ccmplexity 
of the computer model. After the projectile stopped, the 
transiticr frem the force balance equation cn the projectile 
to an equaticn describing the fluid velocity in resfonse to 
the pressure transients was awkward. No simple differential 
equation describes the complexity of the wave mechanics 
invclved with the reflection of the pressure waves in the 
systen. In this regard several alternate apprcaches 
involving sequential alterations of the gcverning equations 
were tried with varying degrees of success. 
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IV. DESCRIPTION CF EXPERIMENTAL APPARATUS 


SS SS ee ee ce ==> er cS ES a Ge cag ae eS SSS Se a 


In crder that the comparison of the experimental and 
ccaputer generated data can be fully understood, a brief 
descripticn cf the NPS experimental apparatus and the 
ccnduct cf tke associated experiments is included (See Ref. 
1, pages 17-27). 


The test chamber was fabricated frcem a three-inch O.D. 
Lucite cylirder, 18 inches long, bored to a 20 mm inside 
Ciameter and fitted with alunzinun end caps. The chamber was 
leaded with ¢ Erase slug weighing 93 grams which rode cn two 
graphite filled Teflon sealing rings. The kErass slug, which 
Sinulated the projectile, was cycled from the breech end of 
the chamker tc the barrel end and returned to the _ treech 
end, completing one hypothetical firing cycle. 


Because cf the desire to vary the charge to mass ratio, 
a variakle chamber velocity was necessary. To accomplish 
this with one chanter, a volume control retaining rod was 
designed intc the system. This brass rod, bored to allow 
gas to pass its length, was threaded through a plate which 
was attached to the barrel end cap holding the red in the 
chamber. The rod, which has a Teflon disc cn the end, not 
cnly establisted the volume of the test chamber, but 
provided a tuffer stop for the slug at the end of its 
forward mcticn. Another Teflon buffer was affixed to the 
breech end cap to cushion the slug in return motion. 


The LEG simulator started a simulated firing cycle with 


the slug at the breech end of the empty chamber, as shcwn in 
Fag. (3) - 
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This is the ready-to-Lran position. The simulated 
prcepellant, cistilled water, was then intreduced, ramming 
the slug tc the opposite (barrel) end of the chamber as the 
Chamber was filled. This was accomplished by applyinc gas 
Fressure tc a power piston which drove the injector piston. 
The injectcr piston forced the propellant past a flow check 
valve and into the chamber. This placed the slug in the 
ready-tc-fire cr ir-battery position. In an actual gun, the 
prcpellant sculd ke ignited at this time in the cycle. Due 


to laboratory constraints, an expulsion system was used. 


An IVDI displacement transducer was manufactured and 
mounted next tc the injector piston. The LVDT was attached 
tc the ccnnecting rod, between the power piston and the 
injectcr fistcn. The volume of the liquid being flaced 
under fpresstre during each shot was measured by filling the 
system in the ready-to-ram position and then draining it 
intc a oaqraduated beaker. By measuring the displacement of 
the injectcr piston head with the LVDT, the volumetric rate 
of fluid injection during the ram stroke was obtained. 


Two fressure taps were drilled in the test chamber. 
These taps were located as close to the breech end as 
pessible, cne at 20 degrees frcem top center and the other at 
20 degrees fron tottom center in a counterclockwise 
direction, as seen from the breech end of the chamber. Only 
tke bottcr iccation was used to record data. 


A 4-chanrel Hewlett-Packard 3960 Magnetic Tape Recorder 
was used to FM reccrd desired data during system oferation. 
For each ¢ata run signals were recorded from two Kaman 
diaphragms type 1000 psi) pressure transducers, one 
ccnnected to the breech pressure tap on the test chamber, 
and the cther to the gas injection side of the power piston. 
The pressure signals were processed with a Kaman Dici-Vit 
Reagout Unit which also provided a visual (digital) display. 
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Data from tke LVDT displacement transducer was also recorded 
cn the tare. 


A Brush Fecorder (Mark 280) was used to cbtain a visual 
display cf tke recorded data. By transcriking the desired 
Signals cn the Magnetic Tape Recorder at a tape speed cf 15 
feet per minute and playing them back into the Brush 
Recorder at :-3/4 feet per minute, the time scale of the 
cutput was expanded by a factor of four on the Brush 
recordings (viz., from a real-time maximum cf 200 mm/sec to 
a delayed tinze maximum of 800 mmysec.) 
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As menticned previously, experimental data was taken for 
ta@ gas fressures ranging frem 50 to 220 psig. From this 
ccllecticn of data, one run at 140 psig was selected as 
representative of system performance. It was felt that any 
model which sculd adeguately describe system operaticn at 
this intermediate pressure would be valid for the entire 
range of expected IPG driving fressures. 


Figure (4) is a comparison cf computed and recorded data 
for injectcr pistor displacement. As can be seen, the 
agreement is very good. This is not very significant in 
that all ucdels tried, as well as the original analog 
computer mocel, were able to correctly predict injector 
pistcn mcticc. | 


Figure (5) is a comparison of analytical results and 
experimental cata for breech chamber pressure. The computer 
model fcllows the shape of the experimental curve fcr the 
duration cf the tine that the projectile is in moticn (0-30 
MSEC). It cscillates very rapidly but does not fall cff to 
a zerc value at the end of the projectile moticn. The 
freguency of these initial oscillations is apfroximately 1.0 
kh2 whick is within the ‘frequency range of the Kaman 
pressure transducers used in testing. However, the mounting 
of these transducers within a connecting cavity instead of 
flush with the chamber wall could have led to then 
experiencing a reduced, lagging frequency response. It is 
felt that scme, if not all, of the ccmputed pressure 
oscillaticns must exist as evidenced by the close 


ccerrespendence of the first two peaks in Fig. (5). The 
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accuracy cf the experimental fressure reading could also be 
considered as being reduced by the location of the pressure 
transducer ir the test setup. It is possible that the 
monitcring j;ressure tap was too close tc the end cf the 
kreech chamber and may not have sensed the full chamber 
Fressure during the latter part of the prcojectile's travel 
wren the fluid velccity is greatest. It is felt that to 
adeguately cescrite the pressure decrease toward zero, it 
would be necessary to completely account for exact changes 
in system cecmetry. To do this would require using a 
distributed parameter approach featuring variable lengths 
fcr the in“ectcr chamber and breech chamber. 


The fpresstre transients which occur after the prejectile 
stcps are ccaplex interference phenomena which are not fully 
described by the mcdel. However, the peak pressure and the 
Natural freguency of oscillation of the computer model 
ccupare fevcrably with the experimental data. 
Unfortunately, the damping characteristics cf the model do 
not follcw the experimental data well. The system 
parameters wkich determine these values will be discussed in 


the next secticn. 


It shculd be reiterated that the lumped parameter 
approach is ar approximation. In a distributed parameter 
approach, pressure and velocity would vary ccntinuously with 
distance witkin the system as well as with time. In the 
‘lumped fparaneter approach, these changes are assumed to 
occur only at the input and output of a building Elcock; 
therefore, fressure and velccity are considered constant 
within the building block. This will always lead to some 
discrepancies when comparing model data to experimental data 
taken at at a fixed point. 


Figure (6) shows the simultaneous pressure history for 
several areas of the LPG feed systen. The input tram 
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pressure, injectcr Chamber pressure, connecting line 
pressure, ard breech chamber pressure are grafhed as 
furcticns cf time. 


Figures (7) and (8) show the simultaneous velocities at 
several pcirts as predicted by the computer model. The 
irjectcr fistcn velocity, fluid propellant velocities at the 
beginning ard end of the connecting line, and projectile 
velccity are plotted as functions of time. 


The injectcr and projectile displacements as functions 
of time are included as Figure (9). AS can be seen fron 
Figures (7), (8), and (9), the novement of the injector 
piston is puch mcre rapid than that of the frojectile. Lt 
is felt that this velocity difference is very much a 
furction cf system gecmetry and should not be considered as 
a generalization for all LPG feed systems. The NPS 
experimental apparatus being modeled had only small fressure 
drops betweer injector and breech chambers. In addition, 
the projectile slug was an order of magnitude lighter than 
the injectcr piston. In large scale systems, with 
significant pressure drops and very massive frojectiles, it 
is guite fessible that the injector pistcn's full stroke 
will occur significantly before the projectile has seated. 
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A. PARAFETEES AFFECTING PROJECTILE RAM TIME 


Several farameters can be varied to "tune" the model to 
correct ram time. In equations (1) and (9) (the fcrce 
Ealances cn the injector pistcn and projectile) empirical 
constants ere included which acccunt for any static friction 
and back fressure (PDI = injector back pressure, PDS = 
prejectile tkack fressure). In eguations (4) and (6), 
resistance ccefficients are introduced to account for system 
resistances. Finally in equations (2), (5), and (8), 


prepellant effective bulk modulus values appear. 


The value for the injector and projectile total tLack 
pressure is difficult to determine precisely. Reference 1 
cited an experimentally determined value of 24 psi. The 
analcg ccmputer model presented in Ref. 2 used a value of 46 
EFsi for a ram pressure of 140 psi. The present model uses a 
value of 72z4 gsi for PDI and 0.5 psi for PDS, or a total of 
24.5 psi tack pressure, which agrees with the experizental 
value. The value of system back pressure will vary with 
driving fressure and the configuration of each LPG system's 


injectcr and rreech chambers. 


Some tasic research has been conducted to relate the 
value of tke resistance ccefficients to fluid properties. 
Unfortunately these studies which are descriked in Refs. 4 
and 5, dealt with steady state fluid flow in corstant 
diameter lires without flew restrictions, making their 
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results inaprlicable to the present LPG feed system prcblen. 
The value cf ws can ke related to the damping ratic of the 


"water hammer" pressure oscillations. The value of Bue was 
determined ty an iterative process to obtain the best fit to 


experimental data. Figure (10) shows the dependence of 
ram-time on By and at keeping all cther variables 


ccnstant. as expected intuitively, increasing the fluid 


resistance creates a larger system pressure drop, resulting 
in a lcwer fressure exerted on the face cf the projectile 
and hence slcwer ram times. The slope cf the two curves are 
almost identical so that no significant advantage would be 
achieved ty system designs which try to minimize either RV 


Or av to the detriment of the other. It should be noted, 
however, tkat higher values for a do tend tc slow down the 
Lam time ucre tkan high values for aa Since the magnitude 


of BV would frobakly depend mostly on the pressure drcp at 
2 


~ 


the gun valve, which seals the breech chamber, its cGesign 


should te closely watched to ensure rapid ram times. 
Likewise the system designer will have to pay close 
attention to the design of the piping and valves in a large 
scale LPG feed system to achieve optimum performance. 


The value cf the effective bulk modulus was determined 
from the undamped natural frequency of the "water hammer" 
Fressure cscillaticns as will be discussed in the next 
secticn. Ite bulk modulus is a fluid parameter which 
characterizes the spring effect of a liguid. The kulk 
modulus can te substantially lowered by the elasticity of 
the chamber and connecting line walls and the amcunt of 
entrapped gas present in the propellant fluid. As can be 
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seen from Figure (11), the value of the bulk modulus has 
coly a spall effect on fram time. Between bulk mcdulus 
values of 10,000 and 100,000 psi (the expected operating 
region of an operational system) the rap time is almost 
ccrstant, varying less than a quarter of a millisecond. 
This is an encouraging result since minimizing entrapped gas 
is a difficult and costly design constraint. 
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P. PARAMETERS AFFECTING CHAMBER PRESSURE CSCILLATICNS 


As eéxfplained in Section III, after the projectile has 
stcpped tke chamber pressure can ke described by a second 
order differential equation (E£gn. 10). If the system 
described by this differential equation is underdamped, the 
scluticn te the equation is an exponentially damped sinusoid 
with a characteristic damping ratio and undamped natural 
frequency. 


From Egn. (10), the system damping ratio, ¥ iss: 


_A3 Rvp 
Ao 20K 


If this damping ratio is increased, the peak chamber 
pressure as well as the dissipation period of the system 
Eressure cscillations will be decreased. From a desicner's 
EFOint of view, it would be advantageous for both of these 
gGuantities to be as low as possible. Hence, it would seen 
that an increase in system resistance to enhance damping 
would be desirable. However, as shown in the previous 
section,an increase in the system resistance also increases 
the ram time- an undesirable occurrence. Thus, alterations 
in system ccnfiguration which affect resistance coefficient 
values will have to be accomplished carefully tc insure 


optimal system performance. 


From Egn. (1C), the undamped natural frequency cf the 


pressure cscillations in the LFG system is: 


K 


Pll, 


The primary effect of ullage would be to decrease the value 





WN = 
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of effective tulk modulus and consequently tc decrease the 
natural freguency of pressure oscillations. Hence, should 
it be fcund in future tests that the actual liguid 
prcpellant is highly sensitive to fressure oscillations, 
changes in ullage may be effected to decrease systen 
Tinging. Fortunately, as shown in the previcus section, the 
effect of ullage variations on ram time is pinimal. 
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Several important conclusions can be drawn frcm the 
output of the LPG ram feed ccmputer nodel. Thus, it is 
realized that serious design consideration must be given to 
the values cf fluid fesistances in the feed system in order 
tc minimize raw time and peak pressures, and to optimize the 
damping cf systen pressure transients. Fortuantely, 
ccncommitart mcdel results also show that the effect of 
entrapped gas is sinimal for the overall operation of a LPG 
feed system. The model was unable to account fcr any 
cavitation effects so this is an area that must be a subject 
cf further studies. 


The fact that large pressure transients cccur during the 
prcejectile ramming cycle indicates the difficulties that 
Bight arise with proposed designs for liguid propellant guns 
with variable propellant volumes and projectile 
displacements. The pressure transients experienced in these 
systems will undcuktedly be very ccmplex in nature. 


If furtker improwements are desired in the modelirg of 
LPG feed systems, the next step should be the inclusicn of 
thermodyrapic effects. For any real fluid, resultant feed 
system pressure and velocity changes will cause temperature 
changes which may substantially affect the propellant fluid 
density. Ccnseguently, propellant igniticn characteristics 
at high rates cf fire will undoubtedly be affected. 


The lumped parameter LPG feed system ccmputer mcdel has 


been shown tc te in favorable agreement with experizental 
data. It is felt that this model has been sufficiently 


43 


= 





validiated tc allow its use in more complicated LEG feed 
system designs. It is expected that this model will be used 
by the Naval Crdnance Station, Indian Head, Maryland to 
assist in studying the fluid dynamic characteristics of a 30 
mm scale model LPG feed system becoming operational in June 
19876. It will be necessary to make the obvious changes to 
descrite tke different geometry of the 30 mm scale sodel. 
The additicnal line lengths, and the larger pressure drops 
due to ccmplex valve arrangements in the 30 mm system will 
reguire different values for model back pressure and 
resistance ccefficients. The presence of an accumulator 
near the Eree€ch ckamber will affect the value of the system 
Gamping ratic. It is possible that the different gecmetry 
and mass characteristics of the 30 mm model could result in 
the injectcr piston being slower than the projectile. 
However Ey inccrforating ali of these alteraticns, the 
manner in which each of the system parameters affect total 
system perfcrmance will be able to be predicted Ey the 
versatile lumped parameter computer model developed in this 
study. 


As such, computer simulaticn is a useful design tool. 
However, tc remain useful, it must be viewed in its freper 
Ferspective. Full scale prototype testing is the only 
conclusive method of demcnstrating ordnance system 
performance. Unfortunately, testing is both costly and time 
consuming. Thus, computer simulation, no matter hew simple 
or complex, can ke used to designate critical testing 
instances and to identify those areas where design effcrts 
Will be acst productive. In this way, the time involved in 
optimizing systen performance as well as subseguent 
preducticn ccsts can ke reduced. In explosives research and 
ordnance design, which is still an empirical and necessarily 
hazardous science, computer simulation can be particularly 
useful in specifyirg the correct approach for design testing 
practices. In this respect, such models as the lumped 
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parameter model developed in this study can ke considered as 
critical signgests at some of the many crossroads in systems 
development. However, it must be realized that as such they 
cnly point tke way. The vehicle for arriving at an 
operational system can only be diligent research and 
engineering kased on a measured progression of system 
demcnstraticns overseen by dedicated project managers aware 
of the many sitfalls along the way. 
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APPENDIX A 


COMPUTER PROGRAM LISTING 


The tTasis of the state variable method of systems 
analysis is tke interpretation that the energy state of the 
ccuponents of a system completely describes the condition of 
a system. As defined in Ref. 6, the state of a system is the 
set of variables, the "state variables", which ccntain 
sufficient information about the present ccndition of the 
system tc permit the determination of all future time 
history cf tke system - provided that all future inputs to 
the system are kncwn. Therefore the energy state cf those 
elements which store energy completely describes the systen. 
In the case of the lumped parameter model, these elements 
are the fluic inertia and capacitance and the inertia cf the 
injectcr fpistcn and projectile. The energy state cf these 
“stcrage elements can be described as a function of tine in 
terms of tke state variabies > pressure, velocity, and the 
system infgut pressure. Only those variables required to 
completely specify the state of the system need to be 
included. 


The first step in arriving at the computer program was 
to take eguaticns (1} through (9) and normalize, or solve 
them for the highest derivative. The state variables then 
become the pressure or velocity associated with these first 
derivatives. The computer program was developed by 
ccmtining eguaticns and defining new variabtles as follcws: 

Lon = u, 
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x(2) = 
X(3) = 
X(4) = 
X(5) = 
X(8) = 
¥(9) = 

Eguation ( 


= At mM <2 mM «et mm 
“"-M MDM HW £& A = 


becomes: 
XDOT(1) = =a xc) = SY xc2y + ok Pe - ot Pot 
= -———— —- == mmo come R _-_ —— 
Mp Mp Mp Mp 
In the ccmputer program, PR is represented by X(30) and KFP 


and PDI keéccme program constants, C (7) and C{2Z), 
respectively. The value of C(2) and C(7) are input on data 
cards. Using the fact that eae equation (2) becomes: 

* - K K e 
MDOT) = = XCD - = XG) 


1 
where K is input to the computer program as C(3). By 


ccmbining eguaticns ,(3) and (4) and using the incomrfressible 
flcw assuupticn ae oo a the next equation is derived: 





A 
XDOT(3) = ( X(2) = Rv a X(3) - X(4) ) 
where ay, tecomes Cc (4). Next equation (5) and the 


inccmpressikle flow assumptions ee V and AV=A V 


24 2 5 3%6 
yield: 
A. K A, K 
won = —)— XG) - 2 = Xp 
As Lo A> Lo 


and as tefcre, C(3)=K. Equations (6) and (7) combine with 
the flow velccity eguation a eo e to yield: 


XDOT(5) = 





A 
( XG) =i. Gy ens 
3 2 Ao 


where RY beccmes C(5)- Since ae equation (8) kecomes: 
a 


K 
XDOT(8) = g X(5) - -- X09) 
le L3 
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where K again beccmes C(3). Finally, equation (9) keccmes: 


-K A A 
XDOT(9) = ~~ X(9) + —2 X(8) - —3 Pps 
Ms Ms Ms 


where KFS and EDS become C(8) and C(9) respectively. The 
auxiliary equations: 


XDOT (6) 
XDOT (7) 


X(1) 
X(9) 


are used tc calculate injector piston displacement, X(6), 


and projectile displacement, X(7). 
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